Measurements of time-dependent CP violation in B° — > ip{2S)Kg decays 
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We report improved measurements of time-dependent CP violation parameters for B°(B°) — » 
ip(2S)Kg. This analysis is based on a data sample of 657 x 10 6 BB pairs collected at the T(45) res- 
onance with the Belle detector at the KEKB energy-asymmetric e + e~ collider. We fully reconstruct 
one neutral B meson in the ip(2S)K < g CP-eigenstate decay channel, and the flavor of the accompa- 
nying B meson is identified to be either B° or B° from its decay products. CP violation parameters 
are obtained from the asymmetries in the distributions of the proper-time intervals between the two 
B decays: S^ s)K o = +0.72 ± 0.09(stat) ± 0.03(syst), A^ 2S)K o = +0.04 ± 0.07(stat) ± 0.05(syst). 

These results are in agreement with results from measurements of B° — > J/ipK° . 
PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 



In the standard model, CP violation in B° meson de- 
cays originates from an irreducible complex phase in the 
3x3 Cabibbo-Kobayashi-Maskawa (CKM) mixing ma- 
trix Q. In the decay chain T(4 1 5) -> B°B° -> /cp/tag, 
where one of the B mesons decays at time tcp to a CP 
eigenstate fcp and the other decays at time itag to a fi- 
nal state /tag that distinguishes between B° and B°, the 
decay rate has a time dependence [H given by 



-\At\/r B0 

P(At) = — < 

4r R o 



<1 ■ 



Sf CP sin(AmdAt) 



+ A fop cos(Am d At) j. (1) 

Here Sf CP and Af CP are the CP violation parameters, 
Tgo is the neutral B lifetime, Arrid is the mass dif- 
ference between the two neutral B mass eigenstates, 
At = tcp — tt a g, and the b- flavor charge q equals +1 (—1) 
when the tagging B meson is a B° (B°). For fcp final 
states resulting from a b — > ccs transition, the standard 
model predicts Sf CP — — £/sin20i Q and Af CP ~ 0, 
where = +1 (—1) for CP-even (CP-odd) final states 
and 4>i is one of the three interior angles of the CKM uni- 
tarity triangle, defined as fa = tt - &rg(V* b V t d/V* b V cd ). 
Measurements of CP asymmetries in b — ► ccs transitions 
have been reported by Belle [!, 0] and BaBar [f| . Results 
from our previously published B° — > tjj(2S)Kg analysis 
were based on a 140 fb _1 data sample corresponding to 
152 x 10 6 BB pairs [3] . Here we report new measurements 
with an improved analysis [3] incorporating an additional 
465 ftr 1 data sample for a total of 605 ft)" 1 (657 x 10 6 
BB pairs). 

At the KEKB energy-asymmetric e + e _ (3.5 on 8.0 
GeV) collider Q, the T(4S*) is produced with a Lorentz 
boost of /?7 = 0.425 nearly along the z axis, which is de- 
fined as opposite to the positron beam direction. Since 
the B° and B° are approximately at rest in the T(4S) 



ccnter-of-mass system (cms), At can be determined from 
the displacement in z between the two decay vertices: 
At ~ Azj (Pjc), where c is the speed of light . 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Chercnkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside the coil is instrumented 
to detect mesons and to identify muons (KLM) . The 
detector is described in detail elsewhere @. Two dif- 
ferent inner detector configurations were used. For the 
first sample of 152 x 10 6 BB pairs, a 2.0 cm radius 
bcampipe and a 3-layer silicon vertex detector (SVD-I) 
were used; for the latter 505 x 10 6 BB pairs, a 1.5 cm 
radius bcampipe, a 4-layer silicon detector (SVD-II) , and 
a small-cell inner drift chamber were used [Kj ■ 

We reconstruct ip(2S) mesons in the decay chan- 
nel (I — e or fi) and J/ijjir + Tt~ decay channel. J/ip 
mesons arc reconstructed in the decay channel 

and include the bremsstrahlung photons that arc within 
50 mrad of each of the e + and e~ tracks [denoted as 
e + e~(7)]. The invariant mass of the J jvb candidates is re- 

< 



e+e (7)" 



quired to be within -0.150 GeV/c 2 < M t 
+0.036 GcV/c 2 and -0.060 GcV/c 2 < M M+/i - -m J/v , < 
+0.036 GeV/c 2 , where rrij/^ denotes the world-average 
J/t/j mass and M e + e -^ and M Al + AI - are the re- 

constructed invariant masses of the e + e - (7) and /x + /x _ 
candidates, respectively. For the ip(2S) — > can- 
didates, the same procedure is used. In this case, the 
invariant mass is required to be within —0.150 GeV/c 2 
< M„_+ K - M - m^s) < +0.036 GcV/c 2 and -0.060 

< +0.036 GeV/c 2 where 



'e T e (7) 
„2 



GeV/c 2 < 



V(2S) 



denotes the world-average ip(2S) mass [llj. For 
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the ^(2S) -> J/ipn+n- candidates, AM = M /+i - 7r+T - - 
Mi+i- is required to be within 0.580 GcV/c 2 < AM < 
0.600 GeV/c 2 . To reduce the fraction of incorrectly 
reconstructed ip(2S) signal candidates, we select tt + it~ 
pairs with an invariant mass greater than 400 MeV/c 2 . 
The Kg selection criteria are the same as those de- 
scribed in Ref. [l2j ; the invariant mass of the pion pairs 
is required to satisfy 0.482 GcV/c 2 < M n+n - < 0.514 
GeV/c 2 . 

We combine the i^(2S) and Kg to form a neutral B 
meson. The B candidates are identified using two kine- 
matic variables: the energy difference AE = i^ ms — 
i?boam and the beam-energy-constrained mass M^ c = 

Vi^bc^m) 2 ~ (P™ 18 ) 2 ' where E£™ m is the beam energy in 
the cms, and E c j^ s and p c ^ s are the cms energy and mo- 
mentum, respectively, of the reconstructed B candidate. 
In order to improve the AE resolution, the masses of the 
selected J/ip and ip(2S) candidates are constrained to 
their nominal masses using mass-constrained kinematic 
fits. For the CP asymmetry fit, we select the candidates 
in the AE-M hc signal region defined as \AE\ < 0.03 GeV 
and 5.27 GeV/c 2 < M hc < 5.29 GcV/c 2 . To suppress 
background from e + e~ — > qq (q = u,d,s, or c) contin- 
uum events, we require that the event-shape variable R2 
be less than 0.5, where R2 is the ratio of second to zeroth 
Fox- Wolfram moments (l3| . 

The b flavor of the accompanying B meson is identi- 
fied by a tagging algorithm [14( that categorizes charged 
lcptons, kaons, and A baryons found in the event. The 
algorithm returns two parameters: the 6-flavor charge q, 
and r, which measures the tag quality and varies from 
r = for no flavor discrimination to r = 1 for unam- 
biguous flavor assignment. If r < 0.1, the accompanying 
B meson provides negligible tagging information and we 
set the wrong tag probability to 0.5. Events with r > 0.1 
are divided into six r intervals. 

The vertex position for the fcp decay is reconstructed 
using charged tracks that have a minimum number of 
SVD hits [15J. A constraint on the interaction point is 
also used with the selected tracks; the interaction point 
profile is convolved with the finite S-flight length in the 
plane perpendicular to the z axis. The pions from Kg 
decays are not used for vertexing. The typical vertex 
reconstruction efficiency and z resolution are 95% and 
78 fim, respectively [T2|. The / tag vertex determina- 
tion is obtained with well-reconstructed tracks that are 
not assigned to fcp- The typical vertex reconstruc- 
tion efficienc y a nd z resolution are 93% and 140 /xm, 
respectively [L2\ . After all selection criteria are ap- 
plied, we obtain 1618 and 1202 events for the l + l~ and 
J/ifm + ir~ modes in the AE-M^ fit region defined as 5.2 
GeV/c 2 < M hc < 5.3 GcV/c 2 and -0.1 GeV < AE < 
0.1 GeV, of which 680 and 712, respectively, are in the 
signal region. 

Figure Q] shows the reconstructed variables AE and 
Afbc after flavor tagging and vertex reconstruction. The 
signal yield is obtained from an extended unbinned 
maximum-likelihood fit to the AE-M^ C distribution. We 



model the shape for the signal component using the prod- 
uct of a double Gaussian for AE and a single Gaussian 
for Mbc whereas the combinatorial background is de- 
scribed by the product of a first-order polynomial for 
AE and an ARGUS [l6[ function for M hc . For the 
4>{2S) — > J/i/>7r + 7r~ mode, there is a background com- 
ponent that peaks like the signal (peaking background) 
in the AE-M\ JC signal region. This peaking background is 
mainly due to the J/^i(1270)°, J/V'^*(892)-7r+ and 
J /ipK^TT + Ti~ modes, with no real — > J/i/jir + ir~ in 
the final state. The fraction of such peaking events is es- 
timated to be 1% from the AM sidebands in data. The 
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FIG. 1: (color online), (a) AE distribution within the Mhc 
signal region, (b) Mbc distribution within the AE signal re- 
gion for B° -» ip{2S)K%. The solid curves show the fits to the 
signal plus background distributions, while the dashed curves 
show the background contributions. The small contribution 
from peaking backgrounds is discussed in the text. 



signal and background shapes for each decay mode are 
determined from Monte Carlo (MC) events; these shapes 
are adjusted for small differences between MC and data 
using a control sample of B + — > ijj(2S)K + [l~7| events, 
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TABLE I: Number of signal B candidates (N s i g ) and esti- 
mated signal purity (p) in the signal region after flavor tagging 
and vertex reconstruction. 



Mode 


N sig 


P 


V(2S)(i+I-)Xg 


628 ± 26 


0.92 ±0.01 


ip{2S){J/ip-K + K-)K° s 


656 ± 26 


0.92 ±0.01 



which have a final state similar to the signal but with 
higher statistics. This sample, where no CP asymme- 
try is expected, is also used to check the potential bias 
in the measurements of CP violation parameters. The 
signal yields and purity in the AE-Mb c signal region af- 
ter flavor tagging and vertex reconstruction are listed in 
Table [J We define the purity as the ratio of the signal 
yield to the total number of candidate events in the signal 
region. 

We determine Sf CP and Af CP by performing an un- 
binned maximum-likelihood fit to the observed At distri- 
bution for the candidate events in the signal region. The 
likelihood function is 



£(Sfcp i Afcp ) — U "Pi (^fcp i -4/cp ; Aij), 



(2) 



GeV/c 2 , -0.03 GcV < AE < 0.20 GeV). The PDF for 
peaking background events, 'Ppeak(Ai), is the same as 
■Psig(Ai) with CP parameters fixed to zero. The term 
P i(A<) is a broad Gaussian function that represents an 
outlier component with a small fraction f \. The only free 
parameters in the final fit are Sf CP and Af CP ] these are 
determined by maximizing the likelihood function given 
by Eq.H 

The unbinned maximum-likelihood fit to the 1300 
events in the signal region results in the CP violation 
parameters, 

S i>{2S)K% = +0-72 ± 0.09(stat) ± 0.03(syst), 
Ai>{2S)k% = +0-04 ± 0.07(stat) ± 0.05(syst), 

where the systematic uncertainties listed are described 
below. We define the raw asymmetry in each At bin by 
(N + - N_)/(N + +iV_), where N + (AT_) is the number of 
observed candidates with q = ±1 (—1). Figure [2] shows 
the observed At distributions for q = +1 and q = — 1 
with no requirement on the tagging quality (top), and 
the raw asymmetry for events with good tagging quality 
(r > 0.5) (bottom). 



where the product includes events in the signal region. 
We only use events with vertices that satisfy \At\ < 70 ps 
and £ < 250, where £ is the x 2 of the vertex fit calculated 
only in the z direction. The probability density function 
(PDF) is given by 



Vi = (1 - /ol) J 



/si g Psig(At')i? s ig(At J - At') 

/pcakPpcak(At')i? s ig(At l - At') 

(1 - / sig - jp Cak )P bkg (At')i?bk g (At 
f ol P ol {AU). 



At 1 ) 



d{At') 
(3) 



The signal fraction / s j g and the peaking fraction / pea k 
depend on the r region and are calculated on an event- 
by-event basis as a function of AE and Mb c . The PDF 
for signal events, 'Psig(Ai), is given by Eq.Q]and modified 
to incorporate the effect of incorrect flavor assignment; 
the parameters tpo and Arrid are fixed to their world- 
average values [Tl|. The distribution is then convolved 
with a resolution function i? S j g (At) to take into account 
the finite vertex resolution. The resolution function pa- 
rameters, along with the wrong tag fractions for the six 
r intervals, wi (I = 1,6) and possible differences in wi 
between B° and B° decays (Awi) arc determined us- 
ing a high-statistics control sample of semilcptonic and 
hadronic b — > c decays d, [l2| . The PDF for non-peaking 
background events, 7\kg(At), is modeled as a sum of ex- 
ponential and prompt components and is convolved with 
a sum of two Gaussians, which parameterizes the reso- 
lution function i?bkg(At). Parameters in "Pbkg(At) and 
i?bkg(At) are determined from a fit to the At distribu- 
tion of events in the AE-M^ C data sideband (Mb c < 5.26 




FIG. 2: (color online). The top plot shows the At distribu- 
tions for q — +1 and q = — 1 with no requirement on r. The 
dashed curve is the sum of backgrounds while the solid curves 
are the sum of signal and backgrounds. The bottom plot is 
the raw asymmetry of well-tagged events (r > 0.5, 45% of 
the total). The solid curve shows the result of the unbinned 
maximum-likelihood fit. 

The systematic errors on Sf CP and Af CP , summarized 
in Table [ill arc evaluated by fitting the data with each 
fixed parameter varied by its 1 standard deviation (a) er- 
ror. The MC-dctermined parameters are varied by ±2tr 
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to take into account possible imperfect modeling in the 
MC. We repeat the CP fit procedure with the new value, 
add the differences in S, A quadratically, and then assign 
the result as the systematic error. The largest contribu- 
tion to Sf CP comes from vertex reconstruction (0.026). 
This includes the uncertainties in the interaction point 
profile (the smearing used to account for the B flight 
length is varied by — fO /im and +20 /im), the tag side 
track selection criteria, the helix parameter correction, 
the Ai| range (varied by ±30 ps), the vertex quality 
cut £ (changed to £ < 150 and £ < 500), the Az mea- 
surement, and imperfect SVD alignment. The last two 
are obtained from a study of J/tpKg. Each physics pa- 
rameter [tb°, Am,;) is varied by the error in its world- 
average value [ill ]. Systematic errors due to uncertain- 
ties in wrong tag fractions are estimated by varying the 
parameters wi, Awj in each r region by their ±lcr er- 
rors. Systematic errors due to uncertainties in the reso- 
lution function are estimated by varying each resolution 
parameter obtained from data (MC) by ±lcr (±2<x). The 
AE, Mb c parameters and signal fraction in each r region 
are varied to estimate the systematic errors. No signifi- 
cant bias is seen by fitting a large sample of MC events. 
The systematic errors from uncertainties in the peaking 
background are obtained by varying the peaking fraction, 
shape, as well as its CP asymmetry parameters. The sys- 
tematic errors from uncertainties in the background At 
shape are estimated by varying each background param- 
eter by its statistical error. We also include the effects of 
interference between CKM-favored and CKM-suppressed 
B — > D transitions in the /tag final state [Hj]. We add 
each contribution above in quadrature to obtain the total 
systematic uncertainty. 

We perform various cross-checks for this measurement. 
A fit to the first data sample (SVD-I) results in the CP 
violation parameter, S = 0.97 ± 0.18, which is consistent 
with our previous result [J. A fit to the CP asymmetries 
of the control sample gives the CP violation parameters, 
S = 0.02 ± 0.05 and A = -0.03 ± 0.03, which are consis- 
tent with no CP asymmetry. A fit to the sideband events 
in the B° — ► ip(2S)Kg data sample gives an asymmetry 
consistent with zero (S = 0.02±0.21, A = -0.04±0.10). 
A lifetime fit to B° -» ip(2S)K° s and B+ -► xjj(2S)K+ 
gives t b o = 1.51 ± 0.05 ps and t b + = 1.62 ± 0.03 ps, 



respectively, which are consistent with the world-average 
values. We also examine the CP violation parameters 
separately for the ip(2S) -> l + l~ (S = 0.84 ± 0.13, A = 
0.14 ± 0.09) and ip(2S) J/^tt+tt" (S = 0.61 ± 0.13, 
A = —0.09 ±0.10) decay modes. We find that all results 
are consistent within errors. 

In summary, we have performed improved measure- 
ments of CP violation parameters sin2</>i and Af CP for 
B° -> ip(2S)K^ using 657 x 10 6 BB events. These 
measurements supersede our previous result Q and 
are in agreement with results from measurements of 
B° — > J J_d>K° [lj|. Combining the results from B° — > 
J/ipK° H and B° -»• tp{2S)K° s decays, we obtain a new 
Belle average sin 201 = 0.650 ± 0.029 ± 0.018. 

TABLE II: Systematic uncertainties. 



Parameter 






Vertexing 


0.026 


0.028 


Wrong tag fraction 


0.006 


0.023 


Resolution function 


0.007 


0.005 


Fit bias 


0.012 


0.011 


Physics parameters 


0.001 


0.001 


Peaking background 


0.006 


0.005 


PDF shape and fraction 


0.001 


0.003 


Background At shape 


0.003 


0.003 


Tag side interference 


0.001 


0.036 


Total 


0.031 


0.053 



We thank the KEKB group for excellent operation of 
the accelerator, the KEK cryogenics group for efficient 
solenoid operations, and the KEK computer group and 
the Nil for valuable computing and Super-SINET net- 
work support. We acknowledge support from MEXT 
and JSPS (Japan); ARC and DEST (Australia); NSFC 
(China); DST (India); MOEHRD, KOSEF and KRF 
(Korea); KBN (Poland); MES and RFAAE (Russia); 
ARRS (Slovenia); SNSF (Switzerland); NSC and MOE 
(Taiwan); and DOE (USA). 



[1] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); 
M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 
652 (1973). 

[2] A. B. Carter and A. I. Sanda, Phys. Rev. D 23, 1567 
(1981); I. I. Bigi and A. I. Sanda, Nucl. Phys. B 193, 85 
(1981). 

[3] Another naming convention /3(= <f>i) is also used in the 
literature. 

[4] K. Abe et al. (Belle Collaboration), Phys. Rev. D 71, 
072003 (2005). 

[5] K-F. Chen et al. (Belle Collaboration), Phys. Rev. Lett. 



98, 031802 (2007). 

[6] B. Aubert et al. (BaBar Collaboration), Phys. Rev. Lett. 

99, 171803 (2007). 

[7] This analysis is improved compared to our previous pub- 
lication [J]. The peaking backgrounds are estimated di- 
rectly from the data sidebands. The signal and peaking 
fraction in the time-dependent PDF are determined for 
each r-bin. In addition to the lepton tracks from the J/ip, 
we use the prompt -k + tv~ tracks in the vertex reconstruc- 
tion for the ip(2S)(J/tpTT + TT~)Kg decay mode. Only lep- 
ton tracks from the J/ip were used in the previous anal- 



6 



ysis. 

[8] S. Kurokawa and E. Kikutani, Nucl. Instrum. Methods 
Phys. Res., Sect. A 499, 1 (2003), and other papers in- 
cluded in this volume. 
[9] A. Abashian et al. (Belle Collaboration), Nucl. Instrum. 
Methods Phys. Res., Sect. A 479, 117 (2002). 

[10] Z. Natkaniec (Belle SVD2 Group), Nucl. Instrum. Meth- 
ods Phys. Res., Sect. A 560, 1 (2006). 

[11] W.-M. Yao et al., J. Phys. G 33, 1 (2006). 

[12] K-F. Chen et al. (Belle Collaboration), Phys. Rev. D 72, 
012004 (2005). 

[13] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 
(1978). 



[14] H. Kakuno et al, Nucl. Instrum. Methods Phys. Res., 
Sect. A 533, 516 (2004). 

[15] H. Tajima et al., Nucl. Instrum. Methods Phys. Res., 
Sect. A 533, 370 (2004). 

[16] H. Albrecht et al. (ARGUS Collaboration), Phys. Lett. 
B 241, 278 (1990). 

[17] Throughout this paper, the inclusion of the charge- 
conjugate decay mode is implied unless otherwise stated. 

[18] O. Long, M. Baak, R. N. Cahn and D. Kirkby, Phys. 
Rev. D 68, 034010 (2003). 

[19] See |http://www.slac.stanford. edu/xorg/hfag/ for up- 
dated results. 



